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Abstract 

We apply the so-called conformal subtraction scheme to predict perturba- 
tively exclusive processes beyond leading order. Taking into account evolu- 
tion effects, we study the scheme dependence for the photon-to-pion transition 
form factor and the electromagnetic pion form factor at next-to-leading order 
for different pion distribution amplitudes. Relying on the conformally covari- 
ant operator product expansion and using the known higher order results for 
polarized deep inelastic scattering, we are able to predict perturbative correc- 
tions to the hard-scattering amplitude of the photon-to-pion transition form 
factor beyond next-to-leading order in the conformal scheme restricted to the 
conformal limit of the theory. 
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I. INTRODUCTION 



An important theoretical issue of the perturbative QCD approach to exclusive processes 
is to check its self-consistency by studying perturbative corrections. This task was carried out 
in the modified minimal subtraction MS scheme for the next-to-leading (NLO) corrections 
to the photon-to-pion transition form factor Ref. the electromagnetic pion form factor 

PHlOl, an d the two-photon process 77 — > M + M~(M = tt,K) fll[. The latter was only 
calculated for the special case of an equal momentum sharing distribution amplitude (DA). 
Only in a recent analysis of the electromagnetic pion form factor |10| has the NLO correction 
to the evolution of the DA been completely taken into account. Discrepancies in the one- 
loop approximation of the hard-scattering amplitude for the pion form factor were clarified 
in Refs. The NLO corrections to the pion form factor and to the process 77 — > M + M~ 
are rather large at accessible momentum transfer; however, it turned out that for the equal 
momentum sharing DA the perturbative corrections to the integrated cross section of the 
charged meson production can almost be absorbed into the electromagnetic meson form 
factor (see Ref. Q). 

It was mentioned long ago that conformal symmetry plays an important role for the 



perturbative description of exclusive processes []T2| , p~3|1 . Unfortunately, conformal symmetry 



is broken by the renormalization of the UV-divergencies appearing in the hard-scattering part 
as well as in the DA. There are two different sources for the breaking of conformal symmetry 
induced by: (i) the running coupling and (ii) the renormalization of the DA. While the first 
one is proportional to the /5-function and vanishes for a non-perturbative hypothetical fixed 
point, the latter appears in the minimal subtraction (MS) scheme; fortunately, it is absent 
in a special subtraction scheme, which we will call conformal subtraction (CS) scheme \Ti 



Let us point out that the Brodsky-Lepage-Mackenzie (BLM) scale setting prescription is 
suitable to absorb conformal anomalies, which are in the CS scheme only proportional to /3, 
in the scale setting prescription of the coupling. 

The understanding of conformal symmetry and its breaking in perturbative QCD has 



been applied recently for the NLO calculation of the leading twist-2 Wilson-coefficients |15 



and the flavor singlet evolution kernels JTB| , [T7[ ] relevant for off-forward processes as well as 
to prove the hypotheses of naive non-Abelianzation for the renormalon chain appearing in 
the evolution kernel |T5| . In this paper we pursue the phenomenological consequences of the 
CS scheme for the exclusive processes mentioned above. 

The rest of the paper is organized as follow. In Section |T| we discuss the scheme de- 
pendence of the perturbative corrections and argue that the NLO corrections calculated in 
the MS scheme should be reanalyzed in this special scheme that ensures conformal sym- 
metry. The consequences of the BLM scale fixing prescription and the evolution effects of 
the pseudo-scalar meson DA to NLO are analyzed in the flavor nonsinglet sector in Section 
PT1| . In Sections [TV] and |V] we numerically investigate the NLO corrections for the photon- 
to-pion transition form factor and the elastic form factor, respectively, where we study quite 
different parametrizations of the pion DA. Applying the conformal operator product expan- 
sion, which is valid in the CS scheme, we are able to predict the photon-to-pion transition 
form factor to next-to-next-to-leading order (NNLO) restricted to the conformal limit of 
the theory. At the moment evolution effects can not be taken into account at this order. 
Conclusions are given in Section V|. The decomposition and the conformal partial wave 
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expansion of functions appearing in the hard-scattering part of the electromagnetic form 
factor are listed in two Appendices. 



II. SCHEME DEPENDENCE OF PQCD PREDICTIONS 



In Ref. [R] we introduced the conformally covariant subtraction (CS) scheme, which 
ensures that the underlying conformal structure appearing in exclusive processes is manifest 
in the conformal limit. Calculations of hard-scattering amplitudes and evolution kernels 
are carried out usually in the MS scheme. Both schemes are related to each other by 
a (finite) refactorization. For instance, the photon-to-pion transition form factor F 77r at 
large momentum transfer factorizes in both schemes according to the standard factorization 
scheme (SFA) [^gOf: 



F 77r (u, Q) = T MS {u, x, Q, fi) ® <f) MS {x, fi) 
= T cs (u,x,Q,fi) ®0 cs (x,/i), 



(1) 



dx. 



where the hard-scattering part T(u,x,Q, fi), depending on the kinematical variables u> and 
Q, the momentum fraction x as well as the factorization scale fi, can be calculated in 
perturbation theory. The DA <p(x, fi) satisfy the Brodsky-Lepage (BL) evolution equation 



2 d 



dfi 2 



(x, fi) = V(x, y; a s (fi)) <g> <j>(y, fi). 



(2) 



Of course, the physical quantities are independent of the chosen scheme, which means that 

T cs (oo, x, Q, fi) = T MS (w, y, Q, y) ® B(y, x, fi), (3) 
cs (x, fi) = B-\x, y, fi) ® MS (i/, fx), 

where the B kernels satisfy B{x,z,fi) ® B~ l (z,y,fi) = 8{x — y). The latter transformation 
implies that the evolution kernel of the DA transforms inhomogeneously: 



V cs (x,y) = B-\x,z) 



V MS {z, z') 



B{z',y)- 



d 



dfi 



2 B-\x,z) 



B(z,y). 



(4) 



Perturbative QCD predictions for physical quantities are given as truncated series in 
a s . This necessary truncation induces their scheme dependence. The unknown higher order 
corrections may be minimized by choosing an appropriate scheme and scale. The problem 
of finding such an optimal renormalization scheme can be attacked with the help of the 



extended renormalization group equations introduced by Stuckelberg and Peterman |22 
which is equivalent to previous work given in Refs. 



To our best knowledge, no 
comparable methods were developed to find also the optimal factorization scheme. For a 
given process there may exist physical arguments to favour a special scheme. A further 
requirement should be that the factorization scheme respects the underlying symmetries of 
the theory, in order that no anomalous terms appear either in the hard-scattering part or 
in the BL evolution equation. 
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The more restricted problem to find the optimal scale in a given scheme has been widely 
discussed in the literature and three quite distinct methods have been proposed: the principle 
of fastest apparent convergence (FAC) p3| , |24| , the principle of minimal sensitivity (PMS) 
27]-p0|, and the Brodsky-Lepage-Mackenzie (BLM) []31] scale setting. The application of 



these methods can yield quite different predictions (see for instance the analyses in Ref. 

H). 

Although conformal symmetry holds only true in the hypothetical conformal limit, it 
should be manifested in the maximally possible manner in the full theory. For exclusive 
processes in which only mesons participate such a factorization scheme is, up to the scale 
setting problem, uniquely defined in the conformal limitQ. However, as discussed above, also 
in the CS scheme anomalous terms proportional to the /3-function are left and cannot, as 
we will see below, be uniquely fixed. According to the BLM scale setting prescription, these 
anomalous terms can be absorbed in the scale setting of the coupling, and the perturbative 
series will be formally the same as in the conformal theory. Motivated by our discussion it 
seems worthwhile to reexamine the known NLO corrections to exclusive processes in the CS 
scheme and to employ the BLM prescription for the absorption of the remaining conformal 
anomalies. 



III. EVOLUTION OF THE FLAVOUR NON-SINGLET DISTRIBUTION 

AMPLITUDE 

A. General formalism 

First we discuss the application of the BLM scale fixing prescription in the evolution of 
the DA. For the convenience of the reader we outline the whole formalism for the solution 
of the evolution equation (|2|) in terms of the conformal partial wave expansion. Using the 
MS scheme, the evolution kernel 

V(x, y- a.) = ^W\x, y) + (g) V«(x, y) + ■ ■ • (5) 



was computed perturbatively in one- and two- loop approximation [|34]-|3"7|1 . The one-loop ker- 



nel is diagonal with respect to Gegenbauer polynomials C| (2x — 1) of order k and with index 
3/2. In the two-loop approximation this property is spoiled and because of the complicated 
structure, the moments 



7Jfcn (a a ) = -2/ dx f dyCl(2x - l)V(x, y; a s ) 1 V)y c! (2y - 1), (6) 
Jo Jo iV n 

where the normalization factor is N n = (n + l)(n + 2)/ (4(2n + 3)), cannot be directly calcu- 
lated. Fortunately, one can use conformal constraints to compute the off-diagonal moments 



In the case where baryons are also involved, it is known that the mixing problem in the evolution 
equation cannot be completely solved by conformal constraints [^] , and thus also in this conformal 
subtraction scheme some freedom remains. 
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in a very economical way p8|J16| , p!7| . Indeed, in the MS scheme the off-diagonal moments 
of the matrix 7 are induced by a special conformal anomaly matrix 7 C and the running 



coupling pSyiByiTf : 



a(0+ 7 c (0+2-&(0,7 =0, 
9 



where the matrices a and b have the following elements: 



bkn{l) 



2(k - l)(k + I + 3)5 kn , 

f 2(1 + n + 3)5 kn - 2(2n 




3) if k — n > and even 
otherwise. 



(7) 

(8) 
(9) 



While the off-diagonal terms, induced by the renormalization of the coupling, are predicted 
by conformal constraints (0), the special conformal anomaly matrix 7 C contains new infor- 
mation. As explained in Refs. p9| , |38| . pT5| . pJ| , this anomaly matrix can be calculated directly 
with the help of modified Feynman rules and it reads to leading order (LO): 



7 C (°)(/) = -b(l)f»> + 



(o) 



w. 



where 



W k r, 



A 



kit 



c f 



-4(2n + 3)0 - n)(k + n + 3)x 

A fc „-V>(fc+2)+^(l) 2A kn 

(n+l)(n+2) ^ (k-n)(k+n+3) 





if k — n > 
and even 
otherwise 



(10) 



k + n + 4" 



■0 



k — n 



+ 2ip (k - n) - ij) (k + 2) - ^(1), 



with ip (z) = £ \nV(z) and C F = 4/3. 

Employing the conformal partial wave expansion, which is given in terms of the eigen- 
functions of the LO kernel V^°\x,y), 



<f){x,Q) 



£ - 

fc=o 



x)x 



N, 



Cl{2x-imO kk { l i)\M{P)) 



red 



(11) 



the BL evolution equation (g) can be perturbatively solved to any order |j38| . Note that 
the expectation values of the operators for odd k vanish, which is indicated by the J2' 
symbol. The evolution of the composite operators appearing in Eq. ( |TT1) is governed by the 
renormalization group equation (RGE), which possesses the following triangular form due 
to Poincare invariance: 



d k ~ 2 , 

V-^°kl = -lk{a s {v))O k i - J2 'ikn ( a s(v))Oni. 



(12) 



n=0 



The off-diagonal matrix elements 7^ (with k > n) appear beyond the LO and are scheme 
dependent. In the MS scheme they can be simple obtained from Eqs. (0)-(]TD). 



In Ref. |14j it has been shown that in the conformal limit of the theory the CS scheme en- 
sures the conformal covariance of the renormalized operators and, therefore, their anomalous- 
dimension matrix is diagonal. The transformation to the CS scheme is determined by a 
matrix B, which depends only on j c : 
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1 



B = 

l + JY 

where the operator J is defined by 



JA:- 



1-JY + J{i c JT) - 



if k - n > 



(13) 



2(fc-n)(fc+n+3) 

otherwise. 



(14) 



This transformation cancels the off-diagonal part of 7 MS , however, it induces an off-diagonal 
term proportional to the /3-function: 
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cs 



d_ 

dg 



B- 1 



B. 



(15) 



Hence, in this scheme the off-diagonal part of the anomalous-dimension matrix is propor- 
tional to the /3-function: 



ND _ P_ a 
Ykn ^kni 

9 



at* 



wIutc A,„ = gAg + 0(aj|). 



(16) 



The BLM scale setting prescription [RI] can now be applied to absorb the off-diagonal 
term (O) into the scale dependence of the coupling: 



a s (ii*) = a s (/i) 



A^ln(^ + 
2tt V /i / 



(17) 



Thus, the off-diagonal term in NLO may be expressed as 



7fcn 



2tt 



A)A2 + 



2tt 



2tt 



(o) , 
Tk +■■■ 



(18) 



7fe(a s (MD) -7fc(a«(^)) + 



where jj* kn = yuexp |A^/7^,°' , | is the new scale. To absorb the NNLO corrections in the 
analogous way, it is necessary to introduce a different second scale and so on. After this 
procedure the RGE (|T9| ) in the CS scheme takes the form: 



d k ~ 2 
^°ki = -7fcMM))0« " E ' [7fcK(/4J, as(»Z), 



- 7 fc(a.O*))] OS- 



(19) 



n=0 



All scales /4 n , /i^* , . . . are uniquely determined by the conformal symmetry breaking term 
Afc n . This anomaly arises from two sources: (i) from the renormalization of the coupling 
constant that enters in the anomalous-dimension matrix and (ii) from the renormalization 
prescription in the CS scheme. 

Here one comment is in order. Since the effects of arbitrary renormalization transforma- 
tions proportional to the (3 function will disappear in the conformal limit, the CS scheme 
cannot be fixed in the full theory. Below we will deal with two such CS schemes referred 
as CSi and CSn. In the first one only the off-diagonal part that is related to the special- 
conformal anomaly matrix will be removed by the transformation fllSD, while in the CSn 
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scheme also the off-diagonal terms proportional to (3q appearing in the anomalous dimensions 
are explicitly removed (obviously, the inhomogeneous part in Eq. (15) cannot be avoided, 
and finally terms proportional to /3 appear). Although the CSn scheme induces an addi- 
tional symmetry breaking term proportional to the /^-function in the hard-scattering part, 
which looks artificial, it will be interesting to compare both CS schemes with each other. 

The evolution equation (||]) can be perturbatively solved with the help of the conformal 
spin expansion. Solving the RGE flT^), which is an inhomogeneous partial first order differ- 
ential equation, the solution has been written in a compact form in Ref. and it is valid 
for an arbitrary scheme: 



(x, Q 2 ) = '<fn(x, Q, Q ) expl- 

n=0 { JC < 



Q % n (fi)\ (0\O nn (Q )\M(P)r d . (20) 



The partial waves (p n (%,Q,Qo) now contain perturbative corrections, which are induced 
by the off-diagonal anomalous dimension matrix. They are also known as expansion with 
respect to the Gegenbauer polynomials 

00 (1 — x)x 3 

Vn(x,Q,Qo) = T,- 1 ^CU2x-l)B d -(Q,Q ). (21) 

k=n k 

The matrix B d ™(Q, Qo) diagonalizes the RGE (|l|) of the conformal operators and is given 
by 

B*» = I _L_ = i + £7 ND + C (7 ND £7 ND ) + " " " , (22) 
where the operator L is an integral operator acting on a triangular and off-diagonal matrix: 

£7 fc T = -[ Q %Z D M exp \-l Qd 4 W) - iM ) ■ (23) 
JQo H { A* J 

Here we do not include radiative corrections at the reference point Qo, i-e. B d ^(Q ,Q ) = 
5kn- Notice that in general the matrix B dia is different from the previously introduced matrix 
B. 



B. NLO analysis 



In the MS scheme the evolution of the DA in NLO was analysed in detail in Ref. [[40 
The main feature showing up in this order is the excitation of higher harmonics due to 
the mixing of the operators, which yields logarithmic corrections in the end-point region. 
This logarithmical enhancement is hidden in the expansion of the partial waves, which is, 
corresponding to Eqs. (p2|)-(pT|), given as 

<p n (x, Q,Q ) = ^-^d (2x - 1) + {x,Q,Q ) + ■■■, (24) 



N n ' 2vr 



No 1 _ f a *(Qo) \ Po 

<P$H*,Q,Q„) = £ ■^ F ^C|(2x-i) I -WW T r 

n+2 P + 7n - % 
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The off-diagonal matrix 7 « ND is determined by the conformal anomalies and in the MS 
scheme it can be easily obtained from Eq. (|7|): 

c (°) a u 

7 (DND = / (0) _ (0)X Ikn -Pobkn , , 

Jkn Uk Jn } 2{k-n){k + n + 3)' 1 ' 

Here the matrices b and 7 C are defined in Eqs. (§) and (0), respectively. 

Above we introduced two conformal schemes CSi and CSn which are obtained from the 
MS one by the transformations 

c(0) 

B kn = 5 kn -^— 7fc " for OS* (26) 

2ir 2{k — n) [k + n + 3) 

c (°) a u 

Bkn - 6kn ~2^2(k~n)(k + n + 3) ^ ° Sl1 ' 

respectively. According to Eq. ( |I5|) these transformations imply that the off-diagonal part 
is 7S ND = -A)A<2 [see Eq. ©] with 

a£ = % 7 " ; 7fc " forCS I; (27) 

fcn 2(A;-n)(A; + n + 3) V ; 

a is c (°) 

4£ = /7 w7 7fc " V forCS„. 

fcn 2(A;-n)(fc + n + 3) 

can be absorbed via the BLM scale setting prescription; so that then the partial waves 

read 



N n v ' fc iT +2 N h ^ " k \ a s (Q* kn ) 



i fo-Won 


So 


V a.(0) J 




/W 0) 


(0) 



(28) 

Finally, we study in the MS and in both CS schemes the scale dependence of the quantity 

mQo)= r dx ^M t (29) 

Jo X 

which enters in different exclusive mesonic processes. For convenience we changed the nor- 
malization of the DA, i.e. (f>(x) = f n ip(x)/2yE, where the pion decay constant is f n ~ 0.131 
GeV. Notice that this integral was originally defined to LO; here, however, it also contains 
the higher order effects caused by the evolution. It is very sensitive to the end-point be- 
haviour of the DA and, therefore, it may serve as a measure for the logarithmic corrections 



due to the evolution, which occur in this region [40 



Inserting the conformal spin expansion (^U|) into the integral (|2~9"D provides the represen- 
tation 
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i(Q,Q ) = J2' I n(Q,Qo)exp\ - \ 

n=0 I J Q° 



''^(AiOlO^QoMP))**, (30) 



2(3 + 2n) a s (Q) ~ , 2(3 + 2AQ (1)ND 
g ° } " (n + l)(n + 2) + "IF" ^ (* + !)(* + 2) 7fc » (31) 



1 _ f a s (Qo) \ 
1 V o.(Q) ) 



ft>+7< >-r<°> 



Wo 

■ + ■ 



a (o) (o) 

which will be evaluated numerically by an appropriate truncation of the series. Here the 
reduced matrix elements are the conformal moments of the DA now normalized as: 

(0\O n M\n(P)) red = f 1 dxd(2x-lMx, f i). (32) 

Jo 

It is worthy to mention that in the CS schemes, the direct integration of the RG equation 
(|19|) implies the following form: 

(n + l)(n + 2) V a s (Q n ) J 

where Q n depends on Q and 

_ q.(g) = (n + l)(n + 2) ~ , 2(3 + 2AQ (0) / _ ct.(Q) \ 

a s {Q n )~ 2(3 + 2n) J^ 2 (k + l)(k + 2) lk \ a s {Q% n )J 1 ] 



ft>+7<°>-r«°> 
1 _ f a 8 (Qo) Y 

v 1 V «.«) ) 



00 



a , (0) (0) 



Notice that in the conformal scheme CSn a contribution of the form 1 — a s (Q) / a s (Q) arises 
in the hard-scattering amplitudes, too. Obviously, their expansion gives an a s -suppressed 
term proportional to /3 . 

Taking into account a sufficient large number of terms in the series (|30]) - (|33"D (details 
for the estimate of the accuracy of such an approximation can be found in Ref. [0), the 
relative correction to NLO, i.e. 

1 {QM = 1^Q») ' (35) 

is numerically evaluated, where we consequently took into account the following perturbative 
expansion: 



exp I - I — 7n(/i) 

JQ.n U 



djx , .1 / a s (Q) 



Pa 



Qo yU J \a s (Qo) 



' <*.(Q) _ <*.(Qo) \ (tP _ A_tF 

,2ti 2tt J { p 2(3 fa 



(36) 



Here ft = 11 - 2/3n 7 , ft = 102 - 38/3n/, 7 i 0) = ~C F {3 + l/[(n + l)(n + 2)] - 4^(n + 2) + 
4-0(1)}, and the expression for 7W can be found for instance in Ref. j|T|. Two values of the 
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QCD scale parameter, a rather small one A LO = 100 MeV as well as a rather large one A LO = 
500 MeV are chosen. Since in general non-perturbative (model) calculations are performed 
at a low scale, we choose the reference momentum Q = \/0.5 GeV. In the following it 
is assumed that the evolution equation can be used at this low scale. The evolution runs 
up to the scale Q = \^20 GeV. The number of active quarks is n/ = 3 and we take the 
following set of DAs: asymptotic, two-hump and ip a (x) = T(2a + 2)/T(a + l) 2 [x(l — x)] a 
with a = {10, 1/2, 1/4}, i.e. very narrow and broad ones are also included. The resulting 
relative NLO corrections in percentage are listed for the three considered schemes in Table 









A = 100 MeV 




A = 500 MeV 


DA 


/(Qo.Qo) 


MS 


CSi 


CSn 


MS 


CSi 


CSn 




2.1 


0.7% 


1.0% 


1.0% 


1.6% 


4.2% 


2.2% 


lf &s 


3.0 


-0.1% 


0.4% 


-0.4% 


-1.3% 


2.3% 


-3.0% 


if 1 / 2 


4.0 


-0.8% 


-0.3% 


-2.4% 


-3.1% 


0.7% 


-7.8% 


<^ Z 


5.0 


-1.2% 


-0.4% 


-2.5% 


-5.6% 


-0.6% 


-11.7% 




6.0 


-1.2% 


-2.3% 


-6.6% 


-6.7% 


-1.6% 


-25.% 



TABLE I. The value of the / integral at the input scale Qq = VCK5 GeV and its relative NLO 
corrections at Q = \/20 GeV are listed for different DAs in the MS, CSi, and CSn scheme for two 
choices of the QCD scale parameter A LO = 100, 500 MeV. 



For both the narrow and the asymptotic DA the NLO corrections are rather small in 
all three schemes. The corrections are negative for broader DAs and their absolute value is 
increasing with growing value of 7 LO , however, they remain rather small in the CSi scheme. 
In the CSn scheme the corrections are larger than in the MS scheme, so one may conclude 
that this conformal scheme is disfavoured. Let us mention that in the MS scheme the 
contributions from the off-diagonal part are of the same sign as the a 2 corrections to the 
diagonal part of the anomalous-dimension matrix and that both of them enter in the relative 
NLO contribution with a similar size. In the MS and the CSi scheme, but not in the CSn one, 
a partial cancellation in the off-diagonal terms takes place between the special conformal 



anomaly and the /3q term |[40|| , and that the net contribution is small. Note that for a 
reference momentum of 1 GeV or even larger the size of the NLO corrections decreases. We 
may conclude that in general the NLO corrections due to the evolution of the DA are small, 
except for very broad DA evolved in the CSn scheme. 



IV. PHOTON-TO-PION TRANSITION FORM FACTOR 



From the theoretical point of view the simplest mesonic process is the production of a 
pseudoscalar meson in two-photon collisions 

l*( qi )l*(q 2 )-M(P), (37) 
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since it is purely electromagnetic to LO. As independent kinematical variables we choose 



the negative of the momentum transfer between the photons squared Q 



-q , with q 



(Qi~ 92) /2 and the asymmetry parameter 00 = Pq/Q 2 . In the case that one photon is on mass- 
shell we have \u\ = 1, while for equal photon virtualities \u\ = 0. The dynamical information 
is contained in the amplitude T a p = —ie 2 e a i3^ u qiq2F^M(u,Q 2 ), where the photon-to-meson 
transition form factor F 7 m(w,Q 2 ) is defined in terms of the time ordered product of two 
electromagnetic currents sandwiched between the one-meson state and the vacuum: 



- ie 6ai3nuqiq%F jM (u, Q 



d 4 xe ixg ^M 



(38) 



At large momentum transfer this transition form factor has been measured for ir°, 77, 
and Vj mesons in single antitagged experiments by the CELLO collaboration |53] and, more 



recently, at CLEO f43|,44j|, where the untagged photon is almost real. The photon-to-pion 
transition form factor has been determined for 0.5 < Q 2 < 2.7 GeV 2 [O and 1.5 < Q 2 < 



9 GeV 2 , respectively WM 



For the CLEO data the virtuality of the second photon was 
estimated to be less than 0.001 GeV 2 . The 77 and 7/ transition form factors are known up to 
20 GeV 2 and 30 GeV 2 , respectively g3|JIIl ( the systematic errors become very large with 



increasing photon virtuality). Data at lower momentum transfer are given in Refs. P5| , ^6 |. 

Taking the SFA to LO [ pi)p0"|| , the normalization of the photon-to-pion transition form 
factor is consistent with pion DA's that are not concentrated in the end-point region. Per- 
turbative and non-perturbative corrections to this prediction have been studied in a number 
of papers. Let us only mention that non-perturbative effects have been included in a model 
dependent way by the transverse momentum dependence P?|-[5T1|] or by the sum rule ap- 
proach [pT]-|53[|. All these analyses show that the data can be reproduced by a DA that is 
close to the asymptotic one|: <y? as = 6x(l — x). It is maybe interesting to note that the 
authors in Ref. pO observed that, taking into account the transverse momentum depen- 
dence in the light-cone formalism, quite different model wave functions are consistent with 
the data. In the following we rely on the SFA, which allows us to calculate the perturbative 
corrections in a systematic way fl~9|] . From this factorization procedure it is expected that 
large soft corrections from the transverse momentum dependence in the pion wave function 
would induce large perturbative corrections in the SFA. 

Alternative to the SFA the transition form factor at large momentum transfer can be 
calculated with the help of the operator product expansion (OPE) at light-like distances. 
In a conformal invariant theory, the form of the Wilson coefficients are fixed up to the 
normalization [|53],[)5[]. In Ref. |T3[] it has been shown that this conformally covariant OPE, 
holds true for (3 = in the CS scheme. Taking into account the general structure of this 
conformal prediction, the decomposition of the transition form factor in conformal partial 
waves is given in the conformal limit as: 



2 In the prediction of the transition form factor only the integral Jq dx p(x)/x enters and its value 
extracted from the data is of about 3 or even smaller. Of course, the asymptotic DA provides 3, 
however, we should mention that one cannot conclude in the mathematical sense that its shape is 
convex; even if one takes into account normalization and positivity. 
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~<k 



k=0 
X 2 Fl 



M«.0) = -^E # 5(*+l.* + 2M^))^^ r ^j II ^ TT (39) 



'k + l + \ lk ,k + 2 + \ lk 
2(k + 2 + \ lk ) 



, - t^-J (*(P)\oM\o) red 

2 lk) l+ujj 

The Wilson coefficients c k and the anomalous dimensions jk are known up to order a 2 
from the perturbative corrections to the longitudinal structure function g 1 |5^j57|] . Here the 
reduced matrix elements are normalized according to Eq. ( |3"2"D and are equal to one in 
LO. 

As already mentioned, the pion transition form factor has been measured by single 
antitagged experiments in the region of 0.5 - 9 GeV 2 at CELLO and CLEO, where the 
antitagged photon is almost on-shell, so that we can set u> = 1. Hence, employing the 
relation 2 F 1 (a,6, c|l) = r(c-i)r"(c-b) l ea ds to a considerable simplification of the prediction: 

F„ (u , = i,Q») = g#f ' , r(* + i)r(* + 2)r(2t + 4 + 7t) 



A. NLO corrections 

It is again worthwhile to mention that in the favoured CSi scheme the prediction (flUD 
is exact to NLO, while in the MS scheme and the CSn scheme it is spoiled by the special 
conformal anomaly and a term proportional to (3q, respectively To evaluate the formula 
([401) to NLO in the CSj scheme, the right-hand side (RHS) of Eq. (|40|) is consequently 
expanded up to the order a s . The NLO results in the two other schemes are easily obtained 
by employing the transformations (p6|) . Table || contains the absolute and relative a s 
corrections to the hard-scattering amplitude, where the virtuality of the space-like photon 
is —q\ = 2 GeV 2 and the QCD scale parameter is set to A LO = 0.22 GeV. For simplicity, 
we use the "natural" scale setting prescription fi 2 = —qf in Eq. (EH) [see the \i dependence 
of the coefficient function] and do not discuss the scale setting dependence. The set of DAs 
is the same as before, however, now at the input scale of \/2 GeV instead of 0.5 GeV, with 
the exception that the two-hump function is evolved from its normalization point at 0.5 
GeV to V2 GeV |58|. For the narrow DA the NLO correction is in all schemes surprisingly 



large, namely, between —25 to —31%. For the asymptotic DA the NLO correction is about 
—20% in both the MS and the CSn scheme, while it is only —12% in the CSi scheme. 
A similar reduction of about 10% in the favoured CSi scheme is observed in the case of 
the two-hump DA resulting in negligible small perturbative corrections in this case. For 
convex amplitudes that are strongly concentrated in the end-point region the corrections 
are positive and becoming very large. These results are due to the fact that only the first 
conformal partial wave has a negative contribution. 

Here one comment is in order. For a given DA, the differences of the predictions in Table 
fjdo not reflect the scheme dependence of the NLO result. One has also to take into account 
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DA 


LO 








NLO 










MS 


CSi 


CSn 


MS 


CSi 


CSn 




0.13 


0.09 


0.10 


0.09 


-31% 


-25% 


-28% 




0.19 


0.15 


0.16 


0.14 


-20% 


-12% 


-21% 




0.25 


0.26 


0.27 


0.24 


5% 


11% 


-5% 




0.26 


0.24 


0.26 


0.22 


-10% 


-1% 


-14% 


(p 1 /* 


0.37 


0.78 


0.76 


0.63 


110% 


104% 


70% 



TABLE II. Absolute LO and NLO predictions for the photon-to-pion transition form factor 
and their relative deviation are listed for different DAs. Here the scale Q is set to y/2 GeV and the 
evolution of the DA is taken into account only for the two-hump function. 



the transformation (JTJ) of the DA, which would reduce the observed differences drastically. 
For instance, if we take the asymptotic DA in the CSi scheme, the input in the MS scheme 

m 



IS 



^ MS ( X ) = B <g> 6x(l -x) = 6x(l - x) \ 1 + ^C F 



> I 1 — X\ 7T 2 

In 2 + 2 



x J 3 



(41) 



The transformation from a given scheme to a second one leads to a logarithmic modification 
of the DA in order a s . In fact the problem arises to which scheme the used non-perturbative 
input is related. Although in the sum rule approach radiative corrections are considered 
as unimportant, for this problem in question, it is interesting to study such perturbative 
corrections. It is expected that radiative corrections are minimized in the CSi scheme. 

Now we would like to discuss shortly the radiative corrections in their dependence on uo. 
Since F 17T (uj, Q) is symmetric with respect to uo — > —u, it is sufficient to consider < u < 1. 
From Eq. ([39]) it is obvious that for u < 1 the conformal partial waves for 2 < k are 
suppressed by u k . For u = only the first term contributes. The first moment of the DA 
is uniquely given in terms of the pion decay constant; moreover, the anomalous dimension 
7o = 0, so the radiative correction, which is also independent of the DA, is determined only 
by the a s corrections to Cq. In the conformal limit of the theory, we may conclude from Eq. 
(|39l) that the to dependence of the perturbative correction to the conformal partial waves is 
rather smooth. 

In Fig. |l](a) the SEA predictions to NLO are compared to the experimental data p2} 444 



for the narrow, asymptotic, and two-hump DA in both the MS and CSi scheme. The 
reference momentum for all DAs is now set to Qq = 0.5 GeV and A LO = 220 MeV. As 
pointed out above the corrections in the conformal subtraction scheme CSi are about 10% 
smaller than in the MS one. Caused by the NLO corrections the prediction of the asymptotic 
DA now agrees better with the data, while the prediction of the narrow DA starts to be 
below the data. In both schemes the two-hump DA remains incompatible with the data. 
However, it can be demonstrated that its prediction can be pushed down for a suitable 
choice of parameters. For instance, if we choose the CSn scheme and A LO = 330 MeV the 
predictions of all three DAs start to be compatible with the data as illustrated in Fig. |l|. This 
is qualitatively the same effect as observed in Ref. [|5(J in the context of the model dependent 
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study of the transverse momentum dependence. Since this scheme transformation does not 
look very natural, we cannot draw any conclusions, however, we may get some motivation 
to study perturbative corrections beyond the NLO. Also we would like to note again that 
as discussed above, the differences in the predictions arise mainly from the fact that we 
"forget" to transform the DA. 
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FIG. 1. The SFA predictions for the pion to photon form factor to NLO are shown for the 
narrow (dashed line), asymptotic (solid line), and two- hump (dotted line) DA, where the reference 
momentum is Qo = 0.5 GeV. In (a) the predictions in the MS and CSi schemes are represented 
as thin and thick lines, respectively, where A LO = 220 MeV. In (b) the CSn scheme is chosen and 
A LO = 330 MeV. 



B. Beyond NLO 



As we can imagine from the results given in the previous subsection, for the phenomenol- 
ogy of the considered transition form factor it is an important task to study the corrections 
beyond the NLO. However, such calculations have not been carried out at present. For- 
tunately, we can rely on the conformal prediction (|3^), valid for (3 = 0, where the Wilson 
coefficients and the anomalous dimensions are known up to order a 2 from the results for 
the longitudinal structure function gi in DIS p6| , |5Tf . At NNLO it is known that the result 
(flOl) is also spoiled in the CSi scheme by a term proportional to the /3-function [|59"l , |l8f . For 
these restricted predictions, symmetry breaking terms proportional to the (3 coefficient are 
neglected. Especially, for the asymptotic DA the NNNLO correction is available from the 
Bjorken sum rule. 

The photon-to-pion transition form factor predicted by the asymptotic DA is given by 
the first term of the conformal expansion ([39]): 



Q 2 F~ /7T (uj, Q 2 



V2U 



3 1 + u 




2lu 



1 +UJ 



c {a s 



(42) 



The coefficient co(a s ) is normalized to 1 at LO. For the case that one photon is almost real, 
i.e. uj = 1, we get 
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Q 2 F(1, Q 2 ) = V2Uco(a.) = 0.185 c (a s ) GeV. (43) 

In the conformal limit the first moment and thus also the asymptotic DA do not evolve. 
The predictive power of the conformally covariant operator product expansion (OPE) tells 
us that the coefficient co(a s ) is the value of the Bjorken sum rule, which is calculated up to 
order [60 , 61] . For three active flavours the numerical result reads^ 



Co (a s ) = 1 - ^ - 3.58333 % - 20.21527 (^) ' + O (at) . (44) 

The higher-loop corrections at a scale of Q 2 = 2 GeV 2 , where A LO is assumed again to be 
220 MeV, reduce the LO prediction to about 17% in NNLO and to about 20% in NNNLO 
(the NLO contribution is 12%). 

From the size of the evolution effects arising in NLO we suspect that the NNLO correc- 
tions remain small and can be neglected in a first step. To obtain the NNLO corrections to 
the hard-scattering part for general DAs, Eq. (fiTf ) is expanded up to order a 2 . In the case of 
the two-hump DA (again evolved to a scale of Q = v2 GeV) the NNLO correction remains 
negligibly small, namely about -2%. For the chosen narrow DA the correction decreases 
from about -24% in NLO to about -30% in NNLO. 

We may conclude that in the conformal limit of the theory the perturbative series looks 
very reasonable for DAs that are not ruled out by the data. It would be very interesting to 
calculate the symmetry breaking effects in NNLO, which can be done by taking into account 
only the rif dependent part of the gluon vacuum polarization, which arise from the quark 
loop. Indeed, this calculation has been carried out recently in the MS scheme for the hard 



scattering part [p9jp!8[ , however, the given representation does not allow to transform this 
result to the CSi scheme in a closed form. Moreover, to be sure that the rif dependent term 
belongs to a conformal anomaly proportional to the (3q term, we have to combine the rif 
part of the hard-scattering amplitude with that of the special conformal anomaly to NLO. 
If one had at hand the desired corrections to the conformal prediction fl3~9|), one would be 
able to fix the scale by the BLM scale setting prescription. 



V. PION FORM FACTOR 

The space-like elastic pion form factor has been extracted for 0.25 < Q 2 < 6 GeV 2 from 
the measured cross section of the process j*p — > n + n. The intermediate off mass-shell pion 
was extrapolated to the pion pole |M],^5|. Since this procedure suffers from large systematic 
errors 



55], a direct measurement of the form factor would be very desirable. Indeed, the 
debate concerning the applicability of the SFA to the electromagnetic pion form factors is 
based on these data, which are doubted in Refs. [|^,|^]. In accordance with the counting 



rule, the data behave as 1/Q 2 for Q 2 > 1 GeV 2 . In the SFA the form factor is predicted to 
LO as: 

FAQ) = 2wf ' C *? (Q) W)\ W=f^. (45) 

6Q Z Jo 1 — x 



3 The a^-correction has been estimated to be negative, too [52,3^ 
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FIG. 2. Predictions for the electromagnetic pion form factor at large momentum transfer Q 
for the asymptotic DA (dashed line) and the two- hump DA (solid line) to LO, where Aqcd = 220 
MeV. Data are taken from Ref. |64l and references therein. 



To fit the existing data shown in Fig. |2], the value of Iq should be about 4-5, so that an 
end-point enhanced DA such as the two-hump function seems to be preferred. 

The general factorization formula for the electromagnetic pion form factor at large mo- 
mentum transfer reads 



F*(Q) = (f>(x, fM F , jj, R ) ® T(x, y, Q, fi F , fi R ) ® <j>{y, fi F , fj, R ), 



(46) 



where the scales are distinguished as the factorization scale [l f and the renormalization 
scale fi R . Note that also the DA depends on the renormalization scale fi R . However, for 
simplicity we will not discuss this dependence in detail and take into account evolution 
effects by the solution of Eq. (|2|), in which both scales are identified. The a s correction to 
the hard-scattering part 



T{x,y,Q' 



xyQ 2 



1 + 



Os(/i 



2tt 



Q,fi F ,fiR) + O (afj 



(47) 



was calculated by different authors in dimensional regularization; however, different renor- 
malizaton and factorization prescriptions were applied @-[I(J. The occurring differences are 
clarified in detail in Refs. |8||| . Indeed, if one takes into account the errors, which are pointed 
, the results given in Refs. fHH can be obtained from those in Refs. ||||-H,[10 



out in Ref. 

The latter results are based on the common renormalization operation and the prescription 
also used in the calculation of the evolution kernel (a detailed discussion about these items 
is given in Ref. M): 



T (i) 



T 



C F T F [ 



x,y 



Q, fi F ) + (3 tP(x, y, Q, fi R ) + (C F - C A /2) T FA (x, y), 



[3 + ln(xy)]lnf^ 

Q 



\n 2 (xy) + 3 \n(xy) 



In x 



Iny 



2(1 -x) 2(1-3/) 
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- In — T- — - ln{xy) + - 



Mr, 



6 



(48) 
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T = Li 2 (l - x) - Li 2 (x) + ln(l - x) In 



y 



i-yj 3 

x + y — 2xy) ln(l — x) + 2xy ln(x) 



;i - x)x 2 + (1 - y)y 2 



(x - y) 2 v 

x [ln(l - x) ln(y) - Li 2 (l - x) + Li 2 (x)] I + {x 



x-y 



y}- 



The appearing In a; and In y terms are responsible for large contributions arising from the 
end-point region. Note that the singular end-point behaviour of T FA (x, y) is — ln(xy), which 
is due to the color structure suppressed by 1/iV 2 = 1/9. The appearing poles in (x — y) 
are actually cancelled by zeros, implying that the corresponding term behaves smoothly for 
x — > y (see Appendix |A]) . 

In contrast to the photon-to-pion transition form factor we do not have conformal pre- 
dictions for the electromagnetic form factor that would enable us to predict higher loop 
corrections without explicit calculation. However, to analyse the NLO correction and to 
compare them in different schemes it is possible to employ the conformal spin expansion: 



F n {Q 



2 _2nf 2 C F a s (fiR)_ g , 2(3 + 2m) 2(3 + 2n) 



3Q 2 

x (n(P)\O mm ^F)\0) red 



njm=0 (m + l)(m + 2)(n + l)(n + 2) 



1 + 



2tt 



(49) 

(o\o n M\n(P)) red , 



where the relative NLO correction is given by the moments: 



T«=4/ dx dy(l-x)CU2x-l)T^(x,y)(l-y)Ci(2y-l). 
Jo Jo 



Using the shorthand notation 



V / /(») 
/ /(x) yel = \ x 



the result reads: 



/(*) 



x ' m 



l ix m.*s^Ack2 X -i), 

o x N m 



Tml — CfT^u + PoT^ n + (C F - C A /2) T, 



FA 
ran J 



(50) 



(51) 



(52) 



T , 



, In X \ rel / In X \ rel 
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X I m 



X 1 n 



'-SH< 



1 / In 2 X \ rel 1 / In 2 X \ rel 



x ' m 2 \ x 1 n 
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In X \ rel 



'lnX\rel 5 /lnX\rel 5 /hlX\rel 14 
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(53) 
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7T 
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(55) 
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where AT^ A and AT^ A are the (relative) conformal moments of the functions AT FA (x) 
and AT FA (x,y) defined in Eq. (|A6|). The singular terms in Eq. (|48f) , given by \r^{x)/x 
for i = 0,1,2, provide a ln*(m)-behaviour and are analytically calculated in Appendix [B|. 
The remaining regular part provides AT FA , AT^ A and AT FA , which are power and 1/N% 
suppressed. Therefore, for all DA this part will become negligibly small for growing m and 
n and thus it is sufficient to calculate only the first few moments. The numerical results are 
listed in Appendix [BL too. 

For the chosen set of DA's it turns out that T F gives a moderate negative contribution, 
while T 13 contains a large positive contribution. There are different suggestions to optimize 
both the factorization and renormalization scale setting. Here we will skip the factorization 
scale setting problem (detailed discussions are given for instance in Refs. [p|, [T0| , |68| ) and deal 



only with the renormalization scale setting prescription. This contribution arises only from 
the renormalization of the coupling and is proportional to (3q. Thus, it can be absorbed in 
a quite natural way by the BLM prescription into the scale of the coupling 



„ f 1 /lnX\rel 1 j\nX\ie\ 5l 

"* -* «- = o«p \- 2 (—L - - 2 (—l + 1) < < 56 > 

with T^ n (Q, Q mn ) vanishing. For the "natural" scale setting prescription fj, F = Q, the form 
factor can be written after summation as 



2nC F f la s {Q 

a s (Q* 
Q 2 



FAQ) = " — ^7^ l |i + ^r«(Q, Q)+0 (a 2 s ) \ I(Q, Q ) 2 , (r,7) 



{l + ^ {1 \Q,Q) + 0(a£)}c*{Q,Q c 



where C n (Q,Qo) — ^C F f^I(Q, Qo) 2 /3 contains all evolution effects, Q* = Qe~ A and the 
/ integral is defined in Eq. (£2^). To obtain the NLO corrections as scheme independent as 
possible, one should take into account the higher order corrections due to the evolution in 
the following manner: 



CAQ, Qo) = Cj°(g, Q ) (1 + ^ • • • + O(aj)) (58) 

and truncate then the series in Eq. ( |57| ) in the first order in a s . Fortunately, for a reference 
momentum square of Q > 2 GeV the size of the NLO contributions arising from the evolution 
allows us to use Eq. (p7|) in practice. In Table |inj we give the results for C v {QoiQo)i r ^ 
and A at Q 2 = 2 GeV for the chosen set of DAs and for the three considered schemes. To 
carry out the calculation more easily, the scale setting has been done after summation. 

Finally, we discuss the phenomenological consequences of the perturbative NLO correc- 
tions. Due to the large value of A the scale Q* appearing in the coupling will be very 
small. For instance, at Q = v2 GeV we find for the asymptotic DA Q* w 0.14 GeV and 
for the two-hump DA Q* w 0.08 GeV in the MS and CSi schemes. At this low scale the 
perturbative treatment of the coupling is not valid anymore. However, it was argued that 
the coupling is frozen at low momentum transfer. A simple parametrization of such a frozen 
coupling arises from the idea that the gluon propagator has an effective mass m g that is 
induced by non-perturbative effects []69|-|71|: 
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DA 






MS 

A 




CSi 

A 




CS n 

A 




0.211 


-8.1 


1.60 


-5.7 


1.60 


-5.7 


1.42 




0.431 


-7.8 


2.33 


-5.2 


2.33 


-5.2 


2.00 




0.767 


-1.22 


3.22 


0.82 


3.22 


0.82 


2.62 




0.848 


-5.9 


2.93 


-3.1 


2.93 


-3.1 


2.46 




1.725 


34.4 


5.10 


29.8 


5.10 


29.8 


3.8 



TABLE III. LO predictions and NLO corrections for the hard-scattering amplitude of the 
electromagnetic pion form factor for different DAs, where the evolution is neglected and the fac- 
torization scale is set to Q. The renormalization scale is fixed by the BLM procedure. 

<*s°(Q) = - ^CTT - ( 59 ) 



In the following we assume that is about 0.5 at very low momentum transfer f72,6? 
Employing Eq. fl5"7|) we obtain from Table |T| the NLO prediction for the electromagnetic 
form factor. In the conformal scheme CSi the prediction is Q 2 F n w 0.12 for the asymptotic 
DA and Q 2 F n « 0.32 for the two-hump DA. These results are compatible with the LO 
predictions shown in Fig. 



VI. CONCLUSIONS 

In this paper we have studied the scheme dependence of the photon-to-pion transition 
form factor and the electromagnetic pion form factor in NLO for two different schemes, 
namely, the popular MS scheme and the so-called conformal scheme. In the case of two- 
particle distribution amplitudes, the latter one is uniquely denned in the conformal limit 
of the theory by the requirement of conformal covariance. Obviously, beyond this limit we 
have to deal with an ambiguity that is proportional to the f3 function. Such terms can 
be naturally included in the scale of the coupling by the BLM scale setting prescription. 



However, as we saw in Subsection |IIIB| in the case of exclusive processes this prescription 
is not sufficient to restore the conformal covariance of the perturbative prediction with 
respect to the original representation. We may expect that for so-called reducible but non- 
decomposable representations conformal covariance holds true. 

Instead of dealing with the usual convolution of hard scattering part and DAs we exploit 
the possibility to work directly with the conformal moments, so that the QCD prediction for 
exclusive processes is given as a sum about such moments. The advantage of such represen- 
tation is that the evolution of the moments is easy to handle and the back transformation 
into the x-space by infinite sums over Gegenbauer polynomials can be avoided. Although 
such sums can be calculated with an appropriate accuracy in a straightforward way, the nu- 
merical cancellation of the oscillations requires some care. For the method, which we used 
here, the exclusive predictions by themselfs are given by infinite sums that do not suffer 
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under oscillations and so they can be calculated without difficulties by taking into account 
a sufficient large number of terms. Especially, for the photon-to-pion transition form factor 
this method allows us to immediately use information available from DIS to give a prediction 
beyond the NLO. 

The NLO corrections in the conformal scheme CSi are in general smaller than in the 
MS scheme. This observation supports the general argument that one should choose a 
scheme in which the underlying symmetries of the theory are preserved in the maximal 
possible manner. In the case of conformal symmetry we cannot fix the scheme uniquely, 
since the renormalization of the coupling causes a conformal anomaly proportional to the 
j3 function. As we demonstrated this remaining freedom can provide us huge differences 
for the NLO corrections. In the case of the photon-to-pion transition form factor we saw 
that due to this freedom the differences between the prediction of end-point narrow and 
concentrated DAs are starting to be compatible with the data. However, one has to be very 
careful with some conclusions. Since we did not transform the DAs, the real question arising 
here is: To which scheme belong the non-perturbative input DAs? It seems to us that this 
problem should be considered with more attention. 

It is also desirable to go beyond the NLO to get more insight in the perturbative cor- 
rections to exclusive processes. At least for the photon-to-pion transition form factor the 
conformal techniques can be extended for the calculation of the NNLO corrections to the 
coefficient function in the full theory. 



In this appendix we extract from T FA (x, y), defined in Eq. fl48|) , that part which vanishes 
in the limit x, y — > 0. Moreover, we derive for these terms an expansion which can be 
conveniently used in the convolution of the hard-scattering amplitude with the DA's. The 
most complicated looking term in T FA (x,y) is proportional to l/(x — y) 2 : 



APPENDIX A: DECOMPOSITION OF T FA 



t 



FA 



1 



(x + y — 2xy) ln(l — x) + 2xy \n(x) + 



(1 - x)x 2 + (1 - y)y 2 



(Al) 



(x - y) 2 



x-y 



x [ln(l — x) ln(y) — Li 2 (l — x) + Li 2 (x)] + {x <->• y}. 



Using the simple algebraic decomposition 



(1 — x)x 2 + (1 — y)y 2 1 — x — y (2 — x — y)xy 
(x-y) 3 x-y (x-y) 3 



(A2) 



and the equality 



Li 2 (l 



x) 



Li 2 (x) — ln(l — x) lnx + 




(A3) 



we immediately obtain for t 
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lFA 



x-y 
ln(xy) 



{1 + (1 — x — y) ln(xy)}m 



1 — x 



(x - yf 



2xy + 



i-yj 

(2 — x — y)xy ^ ( 1 — x 



x-y 



2(l-x-y){Li 2 (x)-Li 2 (y)} 
1 



(x - y Y 



x + y - 2xy) (A4) 
1 — x N 



x In [(1 -x)(l- y)] + 2(2 X _ V)Xy [Li 2 (x) - Li 2 (y)] - (x - y) In . — 

x y yi y 

To extract the desired part, one should take into account the following set of algebraic 
identities: 

fix) - f{y) = _ 1 + x + f( x ) + y + f(y) + y 2 f( x ) ~( x ~ y) x y - x2 f(y) 



x — y 

(l-x-y)[f(x)-f(y)] 
x-y 



x 



y 



(x - y)xy 

+ ±x+(l-x)/(x) + ±y+(l-y)f(y) 
x y 
, (1 - 2x)y 2 f(x) T(x- y)xy - (1 - 2y)x 2 f(y) 



(A5) 



(a; - y)xy 



where in the last equation the upper [lower] sign is used for f(x) = ln(l — x) [Li 2 (x)]. Note 
that the third term on the RHS vanishes for x, y — > 0. Taking also into account the remaining 
part of T FA (x,y), we finally obtain after some simple aglebra the desired representation: 



T^(x,y) = -ln(x) 



ln(l -y)+y 



My) 



ln(l — x) + x 



x 



+ 



7T 2 -7 



AT FA (x) = 



+AT' A (x) + AT (y) + AT' A (x, y), 

ln(l — x) + x ^ ln(x) [(1 — 2x) ln(l — x) + x] ^ ^ (1 ~~ 2x)Li 2 (x) — x 



x 



X 



X 



AT^(x, y) = -2 ln(l - x) ln(l - y) + 6" A (x, y), 



9> A (x,y) 



ln(xy) 
(x - y)xy 
1 

(x - y)xy 



y 2 (l - 2x) ln(l - x) - xy(x - y) - x 2 (l - 2y) ln(l - y) (A6) 



y 2 {2(1 - 2x)Li 2 (x) + ln(l - x)} + xy(x - y) 



-x 2 {2(1 - 2y)Li 2 (y) + ln(l - y)} - xy In ^ — 

1 



1 — x 
2/, 



x 



ln(xy) 

(x-y) 2 

( X + y _ 2xy) In [(1 - x)(l - y)\ + 2(2 ~ - ~ y)xy [Li 2 (x) - Li 2 (y)] 



(2 — x — y)xy . / 1 — x 
2xy + ^ In 

x-y \ l -y, 



(x - y) 2 



x-y 



where AT FA (x) and AT FA (x, y) vanish for x, y — > 0. Especially, for AT FA (x,y) this can be 
seen by a power expansion of 9 FA (x,y), namely, 



3 -ij + f 



— ln(xy) 



2(1 - i - ij + f) 3 + 3i + 2ij - 2j 2 



(2 + i) 



(3 + 



xV 3 . 



(A7) 
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It is now also obvious that 9 FA (x,y) is regular at x — y. In the calculation of the NLO 
corrections it is justified to approximate 9 FA (x,y) by the first few terms of this expansion. 



APPENDIX B: CONFORMAL MOMENTS 



Here we list the needed conformal expansions (with respect to the Gegenbauer polyno- 
mials) of terms appearing in the hard-scattering amplitude of the pion form factor to NLO. 
They are obtained by employing the representation 



x(l — x) 2 



C *(2x - 1) = (-l) k 2{3 + 2k ) 4-rx k+1 (l - x) k+1 



x 



N k ~ KV ~" ~' v ~' (Jfe + l)! dx k 
as well as the definition of the 5-function 



6-1 



(Bl) 



(B2) 



In* x 



x 1 k 



de l Jo 



1 dx x e - l Z(1 „ X) cl(2x - 1)| £=0 



& 2(3 + 2k) r 1 
d? (k + l)\ Jo 



dx 



d k 
dx k 



x 



e-1 



x k+l (l -X 



k+l 

|6=0 



= (-D 



2(3 + 2fc)r(l + e)r(fc + 1 - e) 
r(l-e)r(A; + 3 + e) 



J|6=0 



For % = we immediately find 



1 

x 1 k 



-1) 



k 2(3 + 2/c) 
+ l)(fc + 2)' 



Applying the definition of the digamma function ip(x) = -^lnT(x) provides 



In a; 

x /k 
ln 2 x\ 

X Ik 



2(3 + 2A;) 



pty(l)-^(* + l)-^(fc + 3)], 



(k + l)(k + 2) 

-'I'S 12 *'-*^ 11 "*" 11 ' 

+^'(A; + 1) - ^'(A; + 3) 



(B3) 



(B4) 



(B5) 
(B6) 



For large k we find with ip(k + 1) — ip(l) = In (A; + 1) + 0(1/ k) the following asymptotic 
behaviour 



ln*x 



x 



t = <|) t {[2h(t + l)]' + 0(l/*)}, 



(B7) 



which reflects the singular end-point behaviour for x — > 0. 

Closed formulas for ( ^^-^ )fc and ( ln ^ 1 ~^ +x )k can be derived by employing the expansion 
ln ^~ x ^ = — Y^Lo xl / + !)■ After integration over x the summation can be carried out: 
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/ln(l — x) 



2(3 + 2/c) 



E 



7?! 



x //• (fc + 2)! ^(n-A;)! 
In an analogous way one finds: 

■\n(l-x)+x\ 2(3 + 2k) 



B(n + l,k + 3) = - 



2(3 + 2k) 



(k + l) 2 (k + 2) 2 ' 



E 



nl 



(* + !)! ^ k n + 2(n-k)\ 



2(3 + 2/c) 



B(n + 2,k + 2) 



k + 2,1,2 



(fc + l)(fc + 2) 2 (fc + 3) 3 2 V A; + 3, A; + 4 



(B8) 



(B9) 



Finally, the numerical results that are needed for the calculation of the pion form factor 
to NLO will be listed in Tables [IV| and M. Let us remark that the moments AT£r can be 



m 





2 


4 


6 


8 


10 


12 


14 




/ln(l-aO+aA reJ 
\ x 1 m. 


-0.290 


-0.073 


-0.031 


-0.017 


-0.011 


-0.007 


-0.005 


-0.004 




AT FA 


-2.157 


-0.226 


-0.086 


-0.045 


-0.028 


-0.019 


-0.013 


-0.010 





TABLE IV. Numerical results for the first few moments of ( i£LL_^i±£ \ anc [ AT^ A . The 

\ x I m 

latter is the relative conformal moment of AT FA (x), defined in Eq. |A6|, and appears in Eq. (55) 



rcl 



m 
n 





2 


4 


6 


8 


10 





-0.845 


-0.129 


-0.050 


-0.026 


-0.016 


-0.010 


2 


-0.129 


-0.020 


-0.008 


-0.005 


-0.003 


-0.002 


4 


-0.050 


-0.008 


-0.004 


-0.002 


-0.002 


-0.001 


6 


-0.026 


-0.005 


-0.002 


-0.002 


-0.001 


-0.001 


8 


-0.016 


-0.003 


-0.002 


-0.001 


-0.001 


-0.001 


10 


-0.010 


-0.002 


-0.001 


-0.001 


-0.001 


0.000 



TABLE V. Numerical results for the first few relative conformal moment of AT FA (x, y) defined 



in Eq. |A6| . They appear in Eq. (55) 



easily calculated by employing the representation (|A7|) . 
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